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Abstract 
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1. Introduction 

Nonequilibrium systems of diffusing reactants 
are very common in nature. In chemistry almost 
any chemical reaction is a reaction-diffusion sys- 
tem. In physics, the standard examples are anni- 
hilation of electrons and holes moving in a disor- 
dered media, or vortices and antivortices in type 
two superconductors. Examples from other fields 
include: population dynamics in biology, spread 
of epidemics in health science, and group decision 
dynamics in social science. 

It is customary to denote the various types of 
reactants by capital letters, A, B, C, etc., and 
the rates of the reactions by Greek letters fj., A, a, 
etc. Then, for example, a process in which A and 
B annihilate each other at rate A is represented 

symbolically as: A+B 0. Similarly, a process 
where the reaction of A and B produces C, at rate 
/I, is represented hy A + B C. 

The simplest description of reaction-diffusion 
dynamics employs the densities of the reactants 
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as the basic ingredients of the equations of mo- 
tion. For example, the equation describing the 
binary annihilation reaction 

A + A ^ (1) 

is 

^^DV^nA-^^nl, (2) 

where ua is the density of the reactants, and 
D is their diffusion constant. The first term in 
the above equation represents the diffusive be- 
havior of the particles, while the second term ac- 
counts for the interaction. We shall call these 
kind of equations "mean field equations" for rea- 
sons which will be clarified later on. 

The evolution of many nonequilibrium pro- 
cesses are adequately described by mean field 
equations. One of the most remarkable examples 
is the Belousov-Zhabotinskii reaction where 
a mixture of few chemical reactants produces a 
nonequilibrium process which is periodic in time. 

Facing the success the the mean field theory, it 
is natural to ask whether, indeed, it always gives 
an accurate description of reaction-diffusion sys- 
tems. In fact, it is known that the answer for this 
question is negative. Deviations from the mean 
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field theory appear, usually, in systems of low di- 
mensionality. Returning to the example it 
can be easily seen that the homogeneous solu- 
tion of equation (^) behaves asymptotically as 
ua ~ independent of the dimensionality of 
the system. However, the true asymptotic be- 
havior of (|l|) is UA ~ 1/v^ in Id, tia ^ Int/t in 
2d, and ~ for d > 2 g. 

d = 2 is the critical dimension for reaction- 
diffusion type of nonequilibrium processes. The 
qualitative explanation for this behavior is clear: 
In order to react, the two particles should first 
diffuse to make contact. This diffusion time re- 
stricts the rate of the reaction since, in low di- 
mensions, diffusion is inefficient in mixing the re- 
actants. The Ovchinnikov-Zeldovich segregation 
phenomenon |^ is the result of spatially inactive 
regions developed in diffusion limited reactions. 

The purpose of this work is to present an ex- 
ample for nonequilibrium process in which the 
discretized nature of the reactants has a strong 
impact on the behavior: The mean field equa- 
tions fail to describe the dynamics of the system 
in the long time limit. Discretization effects lead 
to a different asymptotic behavior of the system 
in the long time limit. This change is, in a sense, 
analogous to a quantum phase transition. 

The prototype example we shall use is one of 
the simplest models in population biology: the 
Lotka-Volterra system [^,||. In the predator- 
prey version of this model, two species, a preda- 
tor (A) and a prey (B) are interacting while all 
other environmental factors are assumed intact. 
In the absence of predator, the prey population 
grows exponentially, while in the absence of prey 
the predator death rate results in an exponen- 
tial decay of their population. Binary interaction 
between the species involves the growth of the 
predator population due to consumption of the 
prey, thus: 

A B^2B, and A + B^2A, (3) 

where jU is the predator death rate, a is the prey 
birth rate, and A is the probability for a predator 
to eat a prey at the same spatial location. We 
assume that birth of a new predator follows an 
"eating event" at the same site. The mean field 
equations of this system are the Lotka-Volterra 



equations: 

duA „„2 A 

— — = DV HA - fiUA + Xuaub, 
dt 

= DV'^nB + (TUB - XnAUB- (4) 
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Figure 1. Examples of phase space trajectories 
corresponding to the homogeneous solution of the 
mean field equations (4). 



Here ua and ub denote the population densities 
of the predator and prey, respectively, and it is as- 
sumed that the diffusion constants of both species 
equal to D. The generic behavior of this system 
is periodic in time, see Fig. |l|: When the number 
prey is large, the predator population is grow- 
ing due to the availability of food, but then the 
prey population decreases. Consequently, also 
the predator population diminishes. When the 
predator population is already small, the num- 
ber of prey, again, begins growing and the cycle 
repeats. 

Here, we will show that the behavior of the 
quantized version of this model, in two dimen- 
sions, does not follow the Lotka-Volterra equa- 
tions (^. Our analysis will proceed in the fol- 
lowing way. First, we write down the exact 
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Master equations of the quantized version of the 
system. Then, we map these equations onto a 
a Schrodinger equation in imaginary time, and 
identify the corresponding many-body Hamiho- 
nian. Next, we express the propagator of the sys- 
tem as a field integral and find the corresponding 
action. It will be shown that the mean field equa- 
tions are the saddle point equations associated 
with this action. The effective action of the sys- 
tem will be, then, constructed following the tra- 
ditional procedure of renormalization. Namely, 
the fields will be separated into "fast" and "slow" 
components, and the fast components will be in- 
tegrated out. Finally, we analyze the saddle point 
equation of the effective action and characterize 
the long time asymptotic behavior of the system. 
Our results will be verified by numerical simula- 
tions. 

To avoid cumbersome algebraic manipulations, 
our discussion will be switching between two pro- 
cesses: the binary annihilation (^, and the Lotka- 
Volterra reaction (||). The first example will be 
used as a simple illustration of the derivation, 
while the results for the Lotka-Volterra reaction 
will be usually stated without an explicit deriva- 
tion. 

2. The Master Equations 

The Master equations of a nonequilibrium pro- 
cess are equations for the probabilities of the var- 
ious states of the system. Consider the binary an- 
nihilation process (|^) in zero dimension, i.e. when 
the system consists of a single site. Then, a state 
of the system is defined by the number of reac- 
tants, n, and we denote the probability to find 
the system in this state by P„ . The Master equa- 
tions relate the change of probability in time to 
the rate of flow into and out of the state: 

^ = -| [n(n-I)P„-(n + 2)(n+l)P„+2] . (5) 

The first term on the right hand side represents 
the flow out of the state with n reactants. It 
comes from pair annihilation, and therefore pro- 
portional to the number of pairs, n(n—l)/2. The 
second term accounts for the flow into the state 
which is due to pair annihilation in the state with 
n + 2 reactants. 



The same logic can be used in order to 
construct the Master equations for the Lotka- 
Voltcrra process (||). Considering again the zero 
dimensional case, they take the form 

dPm n 

— -^j^ = —{^J,m + an+Xmn)P„i^n (6) 

+ ld{m+l)Pjn+l,n + <^{n~l)Pm,n-l 

+ A(TO-f)(n + l)P„„ l,n+l, 

where Pm,„ denotes the probability to find the 
system in a state with m predator and n prey. 

The generalization of the above equations to 
the non-zero dimensional case is straightforward. 
A state of the system is now defined by two vec- 
tors of integer numbers: n = (ni,n2,---), and 
m = (rrii, m2, • ■ •). The components of these vec- 
tors represent the occupation numbers of the prey 
and the predator at the various sites of the sys- 
tem, and Pm,n is the joint probability of the oc- 
cupation configurations. The Master equations, 
in this case, contain an additional hopping term 
between the sites. This term will be added to our 
theory later on. 

3. Mapping the Master equations onto a 
Schrodinger equation 

We turn now to map the Master equations onto 
the Schrodinger equation in imaginary time, and 
to identify the corresponding many-body Hamil- 
tonian. Beginning with the example of binary 
annihilation process in zero dimensions, following 
Refs. we define the wave function 

oo 

|V)=^P„|n), (7) 

where \n) denotes a state with n reactants in the 
system. Taking the derivative of \ip) with respect 
to time and substituting (|^) we obtain 

n 

n 

Let us now introduce the creation, a) , and anni- 
hilation, a, operators which satisfy the Bose com- 
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mutation relation [a, a^] = 1, and 

d^|n) = |n+l), while a\n) = n\n — 1) . (9) 

It is easy to see that a^l^j) — J2n ^n+2(" + 2)(n + 
l)|n), and (a^ fa^l^p) = ^„ - l)n|n). Sub- 
stituting these results in equation (|^), one can 
write it in the from of a Schrodinger equation in 
imaginary time 

= (10) 
where the Hamiltonian is given by 
H = ^(atflt _ (11) 

Turning to the Lotka-Volterra reaction in zero 
dimension, we denote by \m,n) the state with m 
predator and n prey. The corresponding wave 
function is \ip) = J2n m Pm,n\'m, n), and the equa- 
tions of motion take the same form as ([l^), but 
with the Hamiltonian: 

H = fi{a^ -l)a+a{l-b)b^b+Xa\b^ -d^)ab. (12) 

Here a) and a are the creation and annihilation 
operators of predator, while 6^ and b are the cre- 
ation and annihilation operators of prey. 

The generalization of equation (|l^) to the 
nonzero dimensional case is obtained by defining 
the creation and annihilation operators at each 
site of the system (i.e. '^Ii'^ ~* Oi,^^ 

bl,b ^ 6i,with [aj,d^] = (5y, etc.), and adding 
a hopping term to the Hamiltonian. The wave 
function, in this case, is 

iv> = E^™.nn(«lr'(^lr'|o)' 

n,m i 

where nii and n,; are the components of m and 
n, and |0) denotes the vacuum state with no re- 
actants in the system. 

4. The formal solution of the Schrodinger 
equation, and expectation values 

The formal solution of the Schrodinger equa- 
tion (^ is 

\m)^uit)\m), 



where ip{0) is the initial wave function, and U{t) 
is the propagator of the system for time t, i.e. 

?7(t) = Tcxpj-^ dt'H{t')^, 

T being the time ordering operator. This solu- 
tion of ip, as function of the time, fully character- 
ize the behavior of the system. Notice, however, 
that the interpretation of the wave function dif- 
fers from that of quantum mechanics. Here it 
represents probability and not probability ampli- 
tude. In particular, the expectation value of an 
operator, Q, in a state defined by ip is given by 
the matrix element 

(Q) - {-PlQm (13) 
where 

i 

is an eigenstate of the annihilation operators, i.e. 

aj\V) ^bj\V = \V) for any j. (14) 

In understanding the structure of this nonequi- 
librium theory, it is instructive to consider spe- 
cific examples of expectation values. Consider, 
first, the expectation value of the identity opera- 
tor Q — 1. Using the identity 

eV(a^a)=/(a^ + l,a)e^ (15) 

where /(d^, d) is a general function of the creation 
and annihilation operators, one can easily see that 
{Vlip) = J2ii m Pm,n- Siucc the sum of probabil- 
ities over all possible occupation configurations 
equals unity, we obtain (Q) = 1. Thus the nor- 
malization of the wave function reads {P\ip) = 1. 

The conservation of probability implies that 
(P|'0) = 1 holds for any time t, therefore 

This equation is satisfied only if {VlH = 0. Thus, 
a legitimate Hamiltonian of our theory must van- 
ish when setting all the creation operators to 
one. For example, the Lotka-Volterra Hamilto- 
nian should satisfy 

i/({d^ = l,6^ = l,d„5J) = 0. 
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It is easy to verify that (|T^) and (|2|), indeed, 
satisfy this condition. 

As a second example, let us calculate the mean 
number of predator at site i, fii = {a\ai). No- 
tice that, unlike quantum mechanics where the 
expectation value of an annihilation operator van- 
ish, here ( [1^ ) implies that fii — {ai). Substi- 
tuting Q ~ hi in ( p^ and using (|l^) we obtain 

~ X^n m ^iPn,m whcre Tii IS the number of prey 
at site i. 

5. The field theoretic formalism 

In order to construct the propagator, U{t), it is 
convenient to employ the path integral formalism. 
To begin with, let us consider the propagator of 
binary annihilation in the zero dimension. We 
define a coherent state of the system as 

|a) ^e-^e-^-^'lO), 

where a is a complex number. The matrix ele- 
ment of a normal ordered operator f{a\ a), where 
all creation operators stand left to annihilation 
operators, is given by 

(a|/(at,a)|a') (16) 

In particular, the inner product of coherent states 
is {a\a') = exp{-i(|ap + |a'p - 2a*a')}, and 
the normalization condition {a\a) — 1 is satis- 
fied. The resolution of identity associated with 
coherent states is 

/ ^ \a){al (17) 

where (Pa = dSia d'^a. 

Consider the expectation value of a general op- 
erator Q = Q{a) , a) at time t, 

(Q(i)) = (P|Q(at,a)e-^(^'''^)'|^), 

where operators are assumed to be normal or- 
dered. For simplicity we choose an initial state, 
1^) = E„^n(a^)"|0), with P„ = e-i/n!, thus 
IV') — e°- ^^|0). Since we are interested in prop- 
erties which are independent of the precise form 
of the initial condition, this particular choice does 
not have an important effect. 



^From ( [l5| ) and general properties of the prop- 
agator it follows that 

m)) = (0|Q(l,a) X 

g-ff(at+l,a)At ^ ^ ^ g-i/(a+ + l,Q)At 

where the square brackets contain a product of 
N infinitesimal propagators for times At = t/N. 
Now, we insert + 1 identity operators ( |l7|) be- 
tween the various terms of the above product. Us- 
ing (|l^) and taking the continuum limit, iV — s- cx), 
we obtain the expectation value {Q{t)) in the 
form of a path integral 

(Q(i)) = J V[a*,a]e'^» g(l, a(t))e-l'^(°)l'+'^* ^ 
where Fq is the action of the system 

Fo^ f dT{a*{T)dra{T)+H[a*{T) + l,a{T)]}, 
Jo 

and I?[a*,a] — Ylri'^'^ ^i''') / ^] ^^'^ measure of 
the integral. 

Having the path integral expression for the 
propagator, the generalization to the finite di- 
mensional case is straightforward. It merely 
amount for the addition of a diffusive term in the 
action. Thus the action of the diffusion-reaction 
process (|^) is 

F^ Jdrdr |a*(a^-DV2)a+|[a*2a2_2a*a2]| , 

where, now, a*(r, r) and a(r, r) are functions of 
the time r as well as the space coordinates r. 
Henceforth we omit the explicit time and space 
dependence of the fields. 

The action associated with the Lotka-Volterra 
process (^ can be derived in the same way, and 
the result is: 

F=J a*{dr-DS/'^)a + b*{dr-D's/^)b+H, (18) 
where 

H = fm*a - (t(1 + b*)b*b + X{a* + 1)(6* - a*)ab. 

Finally we remark that, unlike quantum me- 
chanics of many-body bosons where fields are pe- 
riodic in imaginary time, here there are no such 
boundary conditions. 
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6. The mean field equations 

The mean field equations of tlic field theory 
defined by the action ( p^ ) arc the saddle point 
equations of the bare action, 

SF = 0, 

where the functional derivative is with respect 
to all the fields, a, a*,b and b* . Seeking for a 
solution with non vanishing densities ua — (a), 
fib = (&), the saddle fields (which we denote by 
bar) are a* = b* = 0, and solutions of the equa- 
tions 

da r^„9_ _ , _7 
— = a - + Xab, 
dr 

db o- - - 

— ^ DV^b + ab- Xab. 
dr 

Thus the mean filed equations of the quantized 
Lotka-Volterra system are precisely the Lotka- 
Volterra equations (^. 

The steady state solutions of (|^), duA/dt = 
dng/dt = 0, admit only homogeneous densities 
in space They are associated with the fixed- 
points of the equations: One is the unstable hy- 
perbolic fixed-point, ua = tib = 0, correspond- 
ing to the case with no reactants in the system. 
The second, ua — 0, ub — oo, represents the sit- 
uation where the number of prey grows indefi- 
nitely in the absence of predator. The third fixed- 
point, UA = cr/X, Ub = fJ-/X, is an elliptic fixed- 
point corresponding to a balanced ecological state 
with fixed populations. Linearizing equations 
around the latter fixed-point we obtain: 

l( i:: ) = ( ) . (-) 

where 

Mo-{l 7), (20) 

Aha — riA — tla, and An^ = ub — ti-b- 

The matrix Mq, which we shall call the bare 
mass matrix, is the generator of time evolution 
of homogeneous densities in the vicinity of the 
fixed-point. Its eigenvalues determine the stabil- 
ity properties of the fixed-point. If the real parts 
of the eigenvalues are positive, the fixed-point is 



stable, while if negative, it is unstable. In our 
case, the eigenvalues of Mq are purely imaginary, 
e± = -iziyfjla. It implies that close enough to 
the fixed-point the population densities oscillate 
in time. Moreover, a nonlinear stability analysis 
of the Lotka-Volterra equations shows that the 
system exhibits periodic evolution over the whole 
phase space. The reason is the existence of a con- 
served quantity Q (for homogeneous densities), 

if = -I- - ^ In(ns) - ^ In(n^), (21) 

which confines the phase space trajectories to 
move along concentric closed loops, as shown in 
Fig.0. 

Notice that elliptic fixed-points are unstable 
with respect to small perturbations. Any small 
perturbation might shift the eigenvalues e± off the 
imaginary axis. Consequently the elliptic fixed- 
point will become either a stable or an unsta- 
ble focus. In what follows it will be shown that 
discretization effects, which lie beyond the mean 
field description, indeed, lead to such a scenario. 

7. Renormalization procedure 

An effective action of a field theory is ob- 
tained by integrating out the "fast" degrees of 
freedom. The procedure is usual: First, we 
separate the fields into "fast" and "slow" com- 
ponents. Say for the annihilation process (|^), 
a — ttf + Us, where a/ and Og denote Fourier 
components of the field which oscillate rapidly or 
slowly in space, respectively. Next, we expand 
the action F[a, a*] up to second order in the fast 
fields: F ~ F[as,al]+ F2[as,a*;af,a*j:], where F2 
is quadratic in the fast fields a/ and a^-. Then 
the effective action, Feff[as,a*], is obtained by 
integrating out the fast fields, 

For the binary annihilation reaction, the ef- 
fective action has been calculated by Cardy and 
Tauber |7j. From this calculation they prove that, 
for d > 2, the mean field equations hold, and the 
density of reactants decreases as n(i) ^ 1/t. How- 
ever, when d < 2, inherent spatial fluctuations re- 
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Figure 2. Plot of prey and predator densities as 
a function of time for A = 0.02, ct = 0.006,^ = 
0.008. Upper panel shows simulations with ini- 
tial conditions near the fixed point. Here the sys- 
tem flows away from the fixed point. The lower 
panel shows the flow when the initial conditions 
are set to be far from the flxed point. The con- 
verging behavior of the trajectories signals that 
the dynamics is attracted to a limit cycle. 



suit in a decrease of the asymptotic annihilation 
rate, ~ Xjt'^l'^ 

Notice, however, that in the binary annihila- 
tion example the fixed-point ri^(t ~^ co) — does 
not change due to discretization effects. Only 
the asymptotic approach to this point alters. As 
we shall see, in the Lotka-Volterra system the 
renormalization procedure results in a more dra- 
matic effect. The nature of the fix-point itself 
changes. To put it differently, if one associates 
the long time asymptotic behavior of the system 
with the ground state of the field theory defined 
by the action (jl^), then quantization leads to a 
new ground state which differ from the "classi- 
cal" ground state described by the Lotka-Volterra 
equations (Q). 

In what follows we consider only the physical 
situation of the critical dimension, d — 2, and 
defer the technical details of the derivation of the 
effective action of the Lotka-Volterra system to 
the Appendix. 

Given the fixed-point of the Lotka-Volterra 
equations (jj) at finite densities, {nA,nB) = 
(j,^), it is convenient to change variables to 
fields which represent fluctuations around this 
fixed-point, namely, a a + fiA, and b ^ b + fiB- 
The bare Green function, G°, associated with the 
quadratic part of the resulting action is, 

G" = [{dt - DW^)to + Mor\ (22) 

where tq is the identity matrix, and AIq is the 
bare mass matrix (pO|). Dyson's equation for the 
exact Green function, G, is, 

G = G" + G"I]G, (23) 

where S is the self energy. The correction for 
the mass matrix is the zero Fourier components 
of the self energy, SM = T,{k = 0,uj — 0), thus 
M = Mq + SM. The stability of the flxed point is 
determined by the eigenvalues of the renormalized 
mass matrix, which are 

e± = i (Tr(Af) ± VTr(M)2 - 4Det(M)) . (24) 

Renormalization, in the first approximation, only 
shifts Tr(M) from its zero mean field value. Thus 
the eigenvalues of the mass matrix are approxi- 
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matcly e± « Tr(M)/2 ± iy^JIa, where 

Tr(M) = (6 + log(16)) (a* - a) + (25) 

(tt + 2) ^ + (tt - 2) ^ + (Gtt + 8) 77^5^ I . 

In deriving this resuh (see Appendix) we have as- 
sumed A « i, log (^) » 1, Alog(^) « 

1, and low densities of the reactants. A is the up- 
per momentum cutoff. ^From this result it follows 
that the fixed point becomes an unstable focus. 
This is verified by numerical simulations shown 
in the upper panel of Fig. 

It is natural inquire about the nature of the new 
ground state of the system. What is the global 
behavior of trajectories in phase space when the 
fixed-point is an unstable focus? There are two 
sensible scenarios: (a) One of the other mean field 
fixed-points, ua = ns = or = 0,nB = oo, 
becomes stable due discretization, and all the tra- 
jectories converge to this point, (b) Some other 
type of attractive manifold, such as a limit cycle, 
if formed. In the lower panel of Fig. ^ we present 
numerical results of simulations where we set the 
initial conditions to be far from the fixed point. 
The converging behavior of the trajectories indi- 
cates that the second scenario takes place, and 
the long time asymptotics of the system is that 
of a limit cycle. 

8. Summary 

We have shown that the discreteness of the 
reactants in Lotka-Volterra systems, in two di- 
mensions, results in a behavior (Figs. |^) which 
differs from that of the Lotka-Volterra equations 
(Fig. |l]). Our analytical analysis and numerical 
simulations (Fig. H) indicate that the instability 
of the Lotka-Volterra equations signals the for- 
mation of a new ground state where an attractive 
manifold similar to a limit cycle is formed. Fur- 
ther studies on this subject will be focused on the 
nature of the unstable phase, and the behavior of 
the system in one dimension. 



9. Appendix 

In this appendix we provide some details of the 
summation of the one loop diagrams for the Lotka 
Volterra System. Our discussion will be limited to 
properties of the system near the fixed-point. It is 
therefore convenient to change variables to fields 
which describe fluctuations around this fixed- 
point. Thus changing variables as a ^ a 4- n^, 
6 — > b + fiB, the action (|lj) takes the form 

F = y" drdt a*{dt-D\/^)a+b*{dt-D\/^)b+H, 

where 

H = fib*a-aa*b+a{a*-b*)b*b~aa*a*b 
+ ^i{b* — a*)a* a + X(b* —a*)ba + Xa*ab*b 

— Xa*a*ab+^^{b*a* —a*a* —b*b*). 
A 

To construct the quadratic part of the action 
in the fast variables, i^2, it is convenient to rep- 
resent the fields in terms of their Fourier com- 
ponents, e.g. a(k, w) = / dtdre'^^''^~^^'^*a{r,t). 
Now, we define a vector of fast variables = 
{af,bf,aj,b*^)'^, where the fast fields contain 
terms only with high values of the momentum 
k. Thus: 

where elements of the matrix Z are functions of 
the slow fields (henceforth we drop the subscript 
s of these fields), and: 

= « 

where the bilinear reaction terms have been in- 
corporated into the propagator. 

The integration over the fast variables 
yields (DetZ)-^^^ which can be written as 
exp{— iTr In Z}. Thus the effective action takes 
the form: 

Fef f = ^Tr In G^' + ^Tr ln[l + GqZ] (27) 

The effective mean field equations are obtained 
from the effective action by differentiating with 
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respect to the fields a* and b* . This differentia- 
tion yields the new equations of motion (SF^^f f — 
0): 



6F 1 



Sa* 



Tr 



1 



Go 



6Z 



2 l + GoZ ""Sa 



0, 



=f=0 



SF 1, 



1 



-Tr 

2 1 + GqZ 



G< 



SZ 



0. 



(28) 



a*=b*=0 



Here the trace should be understood as a sum- 
mation both over the indices of the matrices and 
an over A; and uj. The latter summation is loga- 
rithmically divergent, signaling that the long time 
asymptotics of the system is different from the 
mean field equations. In the leading approxima- 
tion we cut off this divergence at the point where 
the perturbation expansion breaks down, namely 
k ~ y\/JiajD. In this approximation the loca- 
tion of the new fixed-point is shifted to: 



rib 



A 
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(29) 



where A is the upper momentum cutoff. Here we 
assume A « l, log ( jg) » 1, A log ( J^) « 

1, and low densities of the reactants, namely 
a/XA^,n/XA^ < 1. 

Next, in order to calculate the mass matrix, we 
expand the equations of motion (|2|) to second 
order in a and b and to first order in X/D. The 
resulting equations are linearized around the new 
fixed point (^ , and the corresponding mass ma- 
trix is calculated. The result for the trace of the 
mass matrix is given in (|26|). 
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